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ABSTRACT 

This  investigation  of  the  radar  cross-section  of  a  finite  cone 
can  be  divided  into  three  areas.  First,  the  exact  solution  for  the 
scattering  of  a  plane  electromagnetic  wave  by  a  finite  cone  is  pre¬ 
sented.  Rigorous  electromagnetic  tneory  is  used  in  the  solution,  and 
no  a]  ocimations  are  made.  Secondly,  methods  of  obtaining  numerical 
results  for  the  radar  cross-section  from  the  analytic  solution  by  using 
a  digital  computer  ore  discussed.  The  third  area  is  a  presentation 
and  discussion  of  the  numerical  results  obtained. 
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THE  THEORETICAL  AND  NUMERICAL  DETERMINATION 

l 

OF  THE  RADAR  CROSS-SECTION  OP  A  FINITE  CONE 

i 

1.  Statement  of  the  Problem 

The  problem  undertaken  la  the  exact  solution  for  the  scattering 
of  a  plane  electromagnetic  wave  by  a  finite  perfectly  conducting  cone. 

We  consider  only  "nose-on"  Incidence  (Figure  l).  In  order  that  the 
entire  surface  of  the  cone  can  be  expressed  as  a  constant  co-ordinate 
surface  In  spherical  co-ordinates,  the  end  cap  of  the  cone  Is  taken  to 
be  e  segment  of  a  spherical  surface  vlth  center  at  the  apex  of  the 
cone.  Time  variations  are  assumed  to  be  given  by  e  and  MJCS  units 
are  used. 

A  considerable  amount  of  effort,  both  theoretical  and  experimental, 
ha 8  been  devoted  to  the  cone  scattering  problem  by  many  workers  No 
attempt  Is  made  here  to  sunmarlze  the  work,  In  view  of  the  very  excellent 
summary  which  appears  In  the  report  by  Kleinman  ind  Senior  (1963)*  It 
should  be  noted  that  the  present  work  Is  an  extension  of  that  done 
earlier  by  Rogers  and  Schultz  (i960},  and  by  Rogers,  Schindler,  and 
Shultz  (1962). 

This  scattering  problem  Is  treated  herein  as  a  boundary-value 
problem  In  electromagnetic  theory  and  no  physical  approximations  are 
mad  .  The  partial  differential  equation  Is,  of  course,  the  vector 
Helmholtz  equation, 

V2  C  ♦  k2  C  -  0  ,  (1) 

.  *4 

where  k  «  2%/\  and  C  may  be  either  the  electric  field  vector  E  or 
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Fig.  1.  Cone  Configuration 


the  magnetic  field  vector  H.  Solutions  of  (l)  are  obtained  in  the 
form  of  infinite  series  containing  unknown  constants.  To  complete  the 
solution  of  the  problem,  these  constants  are  determined  by  satisfying: 
the  necessary  boundary  conditions  for  E  and  H  on  the  surface  of 
the  perfectly  conducting  cone,  the  radiation  condition  at  infinity, 
and  the  finite  energy  condition. 

Numerical  results  have  been  obtained,  and  these  are  compared  with 
experimental  results  obtained  elsewhere,  as  well  as  with  theoretical 
results  obtained  with  the  use  of  approximate  methods. 

2 .  Solution  of  Vector  Helmholtz  Equation 

The  procedure  used  here  for  obtaining  the  solutions  of  the  vector 
Helmholtz  equation  is  well  known  (Stratton,  19^1)* 

Solutions  of  (l)  are 

T  =  7  *  , 

m  «  V  x  (*  r)  ,  (2) 

-10- 
n  -  -  7  x  m  , 

where  r  is  the  radial  vector  in  spherical  co-ordinates  and  4  is  the 
solution  of  the  scalar  Helmholtz  equation 

V2  4  «•  k2  4  =  0  .  (3) 

In  the  region  surrounding  the  cone,  V  .  2  =  v  •  ft  -  o.  Since  V  •  1  +.  o, 

-«  -4 

we  use  only  the  m  and  n  solutions  to  represent  E  and  H. 


4 


It  is  also  veil  known  thnt  the  solution  of  (3)  is 


♦2  (r,«,0)  •  Okr>  0>  [  sin  XS  ]  ’  (U) 

JflV 

o 

where  n  can  have  the  values  1,  2,  3,  or  4  to  represent  Bessel  functions 
of  the  first  kind  (Jn(kr)),  Bessel  functions  of  the  second  kind 
(nn(kr)),  Hankel  functions  of  the  first  kind  (h^(kr)),  and  Hankel 
functions  of  the  second  kind  (h^(kr)),  respectively,  J^J(cos  0)  is 
an  associated  Legendre  function  of  degree  v  and  order  m,  and  we 
let  e  signify  "even"  and  o  signify  "odd"  for  cos  o0  and  sin  o0, 
respectively. 

The  desired  solutions  of  the  vector  Helmholtz  equation  are  then 
obtained  from  (2)  and  (4): 


-n 

m 

e 


ov 


z"(kr)  I^(cos  0)  [ 


sin  o0  "  j 
cos  m0  J  0 


COS  00  "|  - 

sin  n0  J  0 


(5) 


v(  v+1) 

kr 


zv(kr)  0) 


cos  m0  "j  -* 
sin  o0  J  r 


+ 


COS  00  1  -* 

sin  00  J  0 


m 

sin  0 


z^'(kr)  P^(cos 


sin  o0  "]  “• 
cos  o0  J  0 


(6) 


where 


7r  3r  L  r  '"(kr)  ]  ’  “nd  V  V  and  *0  >re  tke 


spherical  unit  vectors. 


3.  Space  Oectlonallzatlon 


s 


One  of  the  most  important  characteristics  of  the  solution  of  this 
problem  is  that  of  dividing  the  space  surrounding  the  cone  into  two 
regions  to  facilitate  the  field  expansions  and  application  of  the 

—4  «4 

boundary  conditions.  The  E  and  H  fields  are  then  expanded  in 

terms  of  the  radial  and  spherical  functions  appropriate  to  each  region. 

Since  the  scattered  fields  must  be  spherically  diverging  waves 

for  large  values  of  the  co-ordinate  r,  the  use  of  llankel  functions  is 

obvious  since  they  possess  the  desired  wave  behavior  as  ~  -  08  .  In 

particular,  since  we  assume  a  time  variation  of  the  fora  eijjt,  the 
U  2 

use  of  z  (kr)  =  h  (kr)  functions  is  necessary  to  achieve  an  outward 
n  nv 

traveling  wave.  At  the  tip  of  the  cone,  however,  the  Hankcl  functions 
possess  a  singularity  the  order  of  which  is  too  large  to  satisfy  the 
finite  energy  condition.  Ibis  characteristic  of  the  radial  functions 
suggests  a  division  of  the  two  regions  at  a  finite  value  of  r. 

The  behavior  of  the  associated  Legendre  functions  indicates  a 
division  of  the  two  regions  at  r  =  b.  This  is  then  the  surface  that 
we  use  to  separate  regions  I  and  II  (Figure  2).  In  region  II,  the 
fields  exist  and  are  bounded  everywhere  in  the  complete  0  domain  of 
0  =  0  to  a,  requiring  the  use  of  only  associated  Legendre  functions 
of  integral  degree.  In  region  I,  however,  0  -  n  is  not  in  the  domain 
of  interest,  allowing  the  use  of  associated  Legendre  iunctions  of  non- 
integral  degree.  It  will  be  seen  that  the  boundary  conditions  will 
determine  the  non-integral  degree  of  each  associated  Legendre  function 


to  be  used  in  region  I. 
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Fig.  2.  Space  Sec tionalization 


The  reader  may  wish  to  refer  to  Rogers  and  Schultz  (i960)  for  0 
more  complete  discussion  of  the  selection  of  the  various  space  divisions 
possible. 


4 .  Field  Expansions 

In  region  I  the  total  fields  are  designated  by  Ej.  and  .  In 
region  II  we  wish  to  keep  the  incident  and  scattered  fields  separate, 

-»i  — »j_ 

and  designate  the  Incident  fields  by  E^  and  and  the  scattered 

fields  by  E^  and  “ir 

In  region  II  the  incident  electric  field  may  be  expressed  (Stratton, 
1941)  as 


ikr  cos  6 
e 


a 

x 


1 

n=l 


o 


n 


-1 

m 

oln 


♦  r 


n 


-1 

n  . 
eln 


(7) 


where 


.  n  2n+l 
Yn  ~  n(  n+l) 


r 

n 


n+l  2n+l 

n(  n+l )  ’ 


(8) 


and  a  is  a  unit  vector  in  the  x  direction.  Ihe  0  variation  in  the 
x 

incident  field  requires  that  m=l  and  forces  us  to  use  odd  m  f  motions 
and  even  n  functions  in  all  expansions  of  the  electric  field. 

The  scattered  field  in  region  II  is  written  as 


\  f  -*4  ,  -*4  ^ 

/'cm,  +  d  n  ,  , 

L,  n  oln  n  eln  y 

n=l 


(9) 


where  c  and  d  are  expansion  coefficients  to  be  determined  frees 
n  n  ^ 

the  boundary  .onditions.  Here  we  have  selected  z^(kr)  =  h^(kr)  and 
the  m  and  n  functions  as  previously  discussed. 


in  region  I  the  total  electric  field  is  expressed  as 

=  )  a  m]L  +  )  b  n\  .  (10) 

I  L  v  olv  U  u  elu 

v  u 

Here  av  and  are  expansion  coefficients  to  be  determined,  and 
u  and  v  are  the  non-integral  degrees  of  the  associated  Legendre 
functions,  which  are  also  yet  to  be  determined. 

The  analogous  representations  for  the  magnetic  field  are  obtained 
from  Maxwell ' s  equations, 

V  x  E  =  -i  go  uq  H  ,  VxH=iGoeQE  ,  (ll) 

and  the  relations, 

V  x  m  =  k  n  ,  ?xn  =  kra.  (12) 


By  using  (7)  through  (10),  in  addition  to  (ll)  and  (12),  and  noting 
that  k  =  oj/u  ,  one  obtains  the  expressions  for  the  magnetic  fields: 


-1 

n  , 
oln 


r 

n 


(13) 


CC' 


n=l 


d 

n 


7 

) 

J 


dM 


)  b  , 

L  u  elu  J 


(15) 


where  T)  is  the  intrinsic  impedance  of  free  space,  vu rJeQ  • 


For  future  reference,  the  field  quantities  are  now  expended  in 


their  entirety: 
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*1  c  [  5  b*  ^(kr)l£(cos  0)  ]  cos  0  or 


+  T  S  a  J  (kr)  PU^C08  ^  +  I  b  J'(kr) 

L  6  0  u  — -  u  u  ^  IT 


cos  0  a 


[  5  aA(kr)  ♦  J  »>  ^(kr)  ^"°°  °_).  ]  sin  0  t  . 

d0  sin  0 


(16) 


^I-  1  {[  rniiT?ii  Jn(Kr)Pi(cos  0) 


cos  0  a 
J  r 


n=  . 


+  [  Y  J  (kr)  Pi(cos  0)  +  r  J' (kr )  ^  ]  cos  0  a 

l  n  n  , "  '■■■„■  n  n  -rr —  j  u 

sin  0  d0 

-  [  Vn(kr)^2  *  rn  Jn(kr>  1  Bl“  0  °0  )  ' 


ff_  =  )  f  r  d  n^V--^  h  (kr)P1(  cos  0)  cos  0  a 

II  L  l  l  n  kr  nx  n'  j  r 


n=l 


+  T  c  h  (kr)  *n^C0S  0)  +  d  h'(kr)  ^^n  "|  cos  0  a 
L  n  n  sin  "6  n  n  do"  J  0 

-  [  =nhn<kr>  ^2  *  d„»n(kr)  ]  aln  0  S  } 

^  *  T[  {  [  §  80  ^Vr  1?  Ju<kr)Ii(cos  0)  ]  Sln  0  ar 

*  [  5  auJ0kr>  -  J  b^(kr)  ^(cos  °>  ]  sin  0  7 


(18) 


+  [  5  auJJ(kr)  Pa^°°9-  •  J  V0kr)  ] cos  f>  a0  } 


(19) 


1^ 


ffII  ■  Tj  1  {  [  Yn  "^ky  ^  JnUr)^(cos  8)  j  =1"  »  ~r 


n=l 


♦  rY  J'(kr)  dIn  -  T  j  (kr)  ^(cos  ®>  1 
L  °  n  d0  n  n  sin  0  ^ 


sin  0  a 


0 


*  [  Y  J’(kr)  ^(cos  8)  .  p  J  ( kr)  din  1  cos  0  s’  }  . 

L  n  n  ^liTe  n  n  ae  J  0  J 


h!t  =  i  ^  ^  c  n^-n-*  ^  h  (kr)P^(cos  0)  1  sin  0  a 
II  i|  L  L  L  n  kr  n  nx  '  J  ^  r 


(20) 


n=l 


♦  c  h'(kr)  -  d  h  (kr)  ^cos  ^  "]  sin  0  a"n 
L  nn  d0~  nn  "iTn  -0-  J  e 

4  cnhn^kr^  Pn^C08  ^  *  d  h  (kr)  ^  ]  cos  0  a\  1  . 

L  n  n  “sirT0  n  n  d0"  J  *  J 

(21) 

Equations  (l6)  through  (2l)  contain  six  sets  of  unknown  constants,  u, 
v,  a^,  b^,  c^,  and  d^.  These  ore  to  be  determined  by  satisfying  the 
boundary  conditions. 


5.  Boundary  Conditions 

We  have  already  satisfied  the  finite  energy  condition  at  the  tip 
of  the  cone  and  the  radiation  condition  at  infinity  by  the  proper  choice 
of  radial  functions  in  each  region.  The  following  boundary  conditions 
remain  to  be  satisfied: 


(a)  [*J] 


r,0 


0  for  0  =  0  ,  r  <  b 

o  — 


<b)  [V^i] 


(22a) 


fra 


for  0  <  0  <  0 


0,0 


0,0 


vr  =  b 


for  0  <  0  <  x 

~  J 


(22b) 
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(c) 

[  %  •  ^  • 

T  Hb  1  for  r  •  b,  0  <  9  <  8 

L  1  J8,0  ‘  ° 

(22c) 

(d) 

The  finite  energy 

condition  at  the  edge  of  the  cone 

(r  -  b,  8  -  eQ), 

( 22d) 

where  b  Is  the  radius  of  the  spherical  cap  and  is 

half  of  the  exterior  apex  angle. 


6.  The  Solution 

To  satisfy  boundary  condition  (22a)  we  first  equate  the  r- component 
of  Ej  to  zero  at  0  =  9q  , 

1  bn  ^Vr~  ^  V'kr)  ^(co‘  8o>  C0B  0  *  0  -  (23) 

u 

and  thus  set 

F^coa  0  )  »  0  .  (24) 

u  o 


This  equation  determines  the  values  of  u.  Equating  the  0- component  of 

to  zero  at  9  =  6  gives 
I  o 


1 


8vJv(‘tr)  — 
v  dd 


l  w(Kr) 


pJ(cos  e0) 
sin  e0 


0  . 


(25) 


Since  P^(co8  0q)  =0  by  (24),  we  set 


4* 


de 


e=e 


(26) 


and  thus  the  values  of  v 


ore  determined. 
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Next,  the  boundary  conditions  (22b)  and  (22c)  are  applied  to 
determine  the  four  sets  of  unknown  expansion  co-efficients  a  ,  b  , 
c  ,  and  d^.  For  the  0  component  of  (22b)  there  results 


GD 

l  C  Vn(kb>  !°4-  -  *  rn  ]  y  0 

sin  e  n  d0  ' 


GO 

+  )  [  c  h  (kb)  ^ cos  ^  ♦  d  h'(kb)  1  co/  0 

L,  1  n  n  '"nl  '  n  n  3T-J  / 
n=l 


^  ,  /bv)  F?;(c08  0) 


a  J  (kb)  W™  P/  ♦  X  \  Jl!(kb)  M  co/0,  0  <  0  <  0  , 
v  V  sin  0  “  U  *  dTJ  /  "  0 


27) 


0  ,  e  <  9  <  *  . 

’  n  —  —  ' 


and  for  the  0  component, 


CD 

•  l  [  Vn(kb)  ^2  ♦  r„  JAtkb>  ^.(°°e  9>  ]  s/^ 

n=l  8in  ® 


CD 

)  I"  c  h  (kb)  ^fn  +  d  h'(kb)  ^n^08  ^  1  si/ 0 

L  d0  n  sin  0  J  7 

'  [  1  avJv(kb)  +  i  »„  j;(k»)  ^(coa  6)  ]  81/0,  0  <  e  <  e  . 
v,  de  sin  0  J  / 


n=l 


r 


(28) 


0  ,  0  <  0  <  x  . 

o  —  — 


Similiarly,  for  the  0- component  of  (22c)  there  results 
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L  Tn',n 


de 


nVkb> 
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00 


n=l 


c  h'(kb) 
n  n 


dF1 

n 


d  h  tut)  pi(cos  o)  i 

n  n  — To  J 


■  J  [  i  •-«<«  £ 

v  d0 


i  w»>  ^  ]  •'/ 


0,  o  <  e  <  9 

—  —  o 

(29) 


and  for  the  0- component. 


00 

i  \  r  ,,/.vv  ^(cos  ©) 
TT  L  L  v;(kb)  -sL — _l 


n=l 


sin  9 


♦  U  [  ‘nbA<*b> 


n=l 


sin  0 


avJ;(kb) 


Py(  cos  0) 
sin  0 


1 


0,  0  <  0  <  0  . 

’  —  —  o 


(30) 


These  four  equations,  (27),  (28),  (29),  and  (30),  are  functions 
of  0,  (27)  and  (28)  over  the  interval  0  <  0  <  n  and  (29)  and  (30)  over 
the  interval  0  <  0  <  '0  .  In  the  solution  of  Rogers  and  Schultz  (i960) 

these  four  equations  vere  manipulated  in  a  process  that  involved 
differentiation  with  respect  to  0.  It  is  well  known  that  an  infinite 
series  can  be  integrated  term-by-term  with  non-stringent  requirements 
on  the  nature  of  convergence,  w.iereas  term-by-tera  differentiation  of 


an  infinite  series  is  valid  only  with  strict  requirements  on  the 
convergence  of  the  series.  Since  the  exact  nature  of  the  convergence 
of  the  infinite  series  expansions  in  (27)  through  (30)  is  unknown, 
we  here  use  an  integration  process,  in  order  to  avoid  the  problems 
encountered  with  differentiation. 


First  we  multiply  (27)  by  F^cos  0),  multiply  (28)  by  sinO^m  , 

m  d0 

and  subtract  the  two  results.  We  then  integrate  the  resulting 

equation  with  respect  to  0  over  the  interval  0  to  n.  It  it  necessary 

to  evaluate  two  integrals  with  limits  of  0  to  it  and  two  integrals 

with  limits  0  to  0q.  The  integrals  are  common  to  boundary  value 

problems  of  this  type  and  can  be  evaluated  by  using  the  associated 

Legendre  differential  equation,  and  (24)  and  (26).  The  integral 

that  appears  as  a  factor  in  the  summation  fortunately  involves 

the  Kronecker  delta,  5  ,  enabling  the  coefficient  c  to  be  separated. 

The  coefficient  d^  is  separated  in  exactly  the  same  manner  except 

that  (27)  is  multiplied  by  sin  0  ^^m  and  (28)  by  P^(  cos  0). 

d0  d/1 

To  separate  a  ,  (29)  is  multiplied  by  sin  9  at  and  (30)  by 

d0 

and  the  results  added.  The  subscript  a  denotes  a  particular 
value  of  the  infinite  set  v.  This  equation  is  then  integrated  with 


respect  to  0  over  the  interval  0  to  0q.  Again  the  integrals  can  be 
evaluated  by  using  the  associated  Legendre  differential  equation,  and 
(24)  and  (26).  Here  the  integral  associated  with  the  av  simulation 
involves  the  Kronecker  delta,  5va  ,  enabling  the  coefficient  aQ  to  be 


separated. 
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Hie  coefficient  b.  is  separated  in  the  same  manner  as  is  a  , 

p  <1  a 

1  dPX 

except  that  (29)  is  multiplied  by  P  and  (30)  by  sin  0  _ £  . 

p  de 

The  subscript  0  denotes  a  particular  value  of  the  infinite  set  u. 

If  the  values  of  and  F^  given  by  (8)  are  then  substituted  in  the 
four  separated  equations,  there  result: 


m 


-1“  (2m*l)  J  (kb) 
_ m _ 

m(m+l)  h  (kb) 
m 


(2m+l)  sin  0 


♦  0  _ 2|  \ 

2rm(nH-l)" 2  h  (kb)  0=0 
L  J  m  o  v 


dP:  >  a  v(vfl)j  (kb)pj;(  COB  0  ) 


0=0  v(v+i)  -  m(m+l) 

o  v  '  ' 


(31) 


m 


lml  ( 2m+l)  J^(kb) 

m(m+l)  h'(kb) 
m 


(2ml)  P  (cos  0  ) 


♦  T - Iz -  i  avJv(kb)Pt(c08  eo> 

2[»(»*1)]2  h;(kb)  ^  V  v  V  0 


(2nr+l)  sin  6  P^cos  8  )  ,  b  J'(kb)  dp1 

— - - 2-  )  (32) 

. ' ( kb)  m(m+l)  -  u(u+l)  0=0 


2m(nKl)  h 


u 


a 


a 


dP1 

“■ !» V  )  { r ,  ,  ■“  la-1)J»(1 “’h 

w««  ,.iu  ■"  ■<■*» 


d0  0=  0O 
a(a*l)-n(  n+l) 


} 


P^( cos  0  )  v-  (  r 

-2 — 2 - s  { r  d  h  (kb)- 

«(«*!)  <W«>> 


l1^1  (2ml) 
n(ml) 


jn(kb)]  ^(cos  8o)  } 


(33) 


sln  ®0  dS^e-e  * 


n-*-l 


be  " 


b(s*D  b  j  <m>)  n=1 


n  \'  r  r  1  n(n+l)P  (cob  9  ) 

f  H  [  W“>  -  1  n  0 

n(  n+ 1 ) 


n(n+l)  -  u(u+l) 

(3»0 


where  the  quantities  IL  and  B0  are  defined  by 

u  p 


r  0  i  2 

bt  =  J  sin  e  (p^)  de  . 


(35) 


The  reader  may  wish  to  refer  to  Appendix  A  for  the  analytic  details  of 
the  derivation  of  (3l)  thru  (  35 ) * 

Equations  (31)  through  (3^)  could  be  manipulated  into  four  equations 
with  each  set  of  coefficients  appearing  in  only  one  equation,  but  the 
form  of  the  end  result  would  be  less  convenient  for  numerical  computation. 
Therefore,  (31)  through  (3*0,  together  with  ( 16)  through  (2l),  (24),  (26), 
and  (35)  represent  the  formal  solution  of  the  problem. 

We  have  completed  the  solution  without  the  necessity  of  satisfying 
boundary  condition  (22d),  the  finite  energy  condition  at  the  edge  of  the 
cone.  Rogers  and  Schultz  (i960)  used  numerical  results  to  show  that  this 
finite  energy  condition  appears  to  be  satisfied  at  the  edge  of  the  cone. 

One  of  the  primary  objectives  of  the  solution  of  this  problem  is 
to  investigate  the  radar  cross-section  of  the  cone.  The  radar  cross- 
section,  0,  is  defined  to  be 


a 


=  lim  4xr^ 

j*  ao 


(36) 


-  1  r  -  -»*  I 

where  S  =  tj  Re  |  E  x  H  j,  the  average  Poynting  vector.  For  our  co¬ 
ordinate  system,  the  radar  cross-section  evaluated  at  9  =  0  is  more 
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precisely  termed  the  back  scattering  radar  cross-section,  o„.  By 

Bo 

using  some  simple  algebra,  o ^  can  be  shown  to  be  expressed  by 

2  |  2 

aBS  =  h  I  1  in  n(n+l)  (cn'idn)  '  {37) 

n 

where  X.  is  the  wavelength  of  the  incident  plane  wave.  In  order  to 

determine  the  back-scattering  radar  cross-section,  then,  we  mast 

first  determine  the  sets  of  c  and  d  . 

n  n 


7.  The  Numerical  Solution 

Equations  ( 31)  through  (3*0  represent  an  infinite  number  of 
equations  in  an  infinite  number  of  unknown  expansion  coefficients. 
The  expansion  coefficients,  therefore,  do  not  enjoy  the  property  of 
finality.  It  is  important,  then,  to  calculate  as  many  of  the  co¬ 
efficients  as  possible  in  order  to  insure  that  the  values  of  the 
lowest  order  coefficients  are  reasonably  accurate.  Tbe  number  of 
coefficients  calculated  in  each  set  is  designated  by  nQ.  All 
numerical  work  was  done  for  =  165°  (o  cone  apex  angle  of  30°). 
The  calculations  have  been  carried  out  for  a  rather  large  number  of 
values  of  ka  in  order  to  determine  rather  well  the  details  of  the 


graph  of  a  _  vs.  ka,  a  being  the  radius  of  the  base  of  the  cone. 

Bb 

An  examination  of  (31)  through  (3*0  indicates  that  the  following 


sets  of  constants  need  to  be  determined: 


u,  v,  B^,  Bv,  P*(cos  165°), 
>  Jn(kb),  J^(kb),  Jv(kb), 
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J'(kb),  J'(kb),  J'(kb),  h  (kb),  and  h'(kb).  The  first  thirty  values, 
n  wb  •'  n  n 

each,  of  u  and  v,  as  determined  from  (24)  and  (26),  were  taken  from 
Waterman’s  paper  (1963) ,  and  these  have  seven- place  precision.  With 
the  exception  of  the  radial  functions,  the  remaining  sets  of  constants 
were  calculated  by  Schultz,  Bolle,  and  Schindler  (1963),  using  a 
Burroughs  Datatron  205  computer.  The  reader  may  wish  to  refer  to 
their  work  for  a  detailed  presentation  of  the  methods  used  in  calcu¬ 
lating  these  constants.  Appendix  B  herein  li6ts  the  first  fifty-one 
values  in  each  set  of  these  constants  along  with  the  given  sets  u  and 
v.  The  radial  function  constants  were  calculated  by  using  the  infinite 
series  representations  for  the  spherical  Bessel  functions. 

When  these  constants  are  substituted  into  (31)  through  (34),  4nQ 
equations  in  4n^  unknowns  result.  These  4nQ  equations  must  then  be 
solved  for  the  four  sets  of  nQ  expansion  coefficients.  Two  different 
methods  were  used  to  accomplish  this.  First,  an  iterative  method  was 
used  for  values  of  ka  in  the  Rayluigh  region.  Secondly,  a  more 
complicated  method,  but  one  that  is  usable  for  any  value  of  ka  (a 
standard  matrix  solution),  was  used  for  the  higher  values  of  ka. 

For  the  iterative  solution  we  first  assume  initial  values  of  the 

c  and  d  in  (33)  and  (34)  and  obtain  initial  values  of  the  a.  and  b  . 
n  n  v  u 

These  values  of  av  and  b^  are  then  substituted  into  (31)  end  (32)  and 

new  values  for  the  c  and  d  are  obtained.  The  process  is  then 

n  n 

repeated,  and  continued  until  the  values  of  the  coefficients  approach 
a  final  value.  This  method  did  not  converge  for  values  of  ka  greater 
than  O.518. 
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The  second  method  is  more  complicated  bat  mere  useful  since  it 
is  applicable  for  higher  values  of  ka .  This  method  involves  a 
straightforward  matrix  multiplication.  If  nQ  coefficients  are  to  be 
calculated,  (3l)  through  (34)  can  be  written  in  the  form 


c  = 


E  +  E0  a 

1  L* 


(38) 


d  =  Fx  F2  a  4  F3  b 


(39) 


a  =  Gi  c  +  G2  +  G3  d  *  G4 


(^0) 


b  =  H1  d  +  H2  , 

where 


*  ■ 

■  « 

C, 

d 

1 

L 

C2 

P-i 

II 

• 

* 

• 

• 

• 

• 

c 

d 

n 

n 

0 

0 

m  * 

•  m 

(41) 


(42) 


»  m 

a 

’b 

V, 

1 

1 

a 

b 

V 

2 

2 

• 

r  = 

• 

a 

t 

V 

u 

n 

n 

0 

0 

m  m 

to  m 

(43) 
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El’  Fl'  G2'  G4'  and  H2  are  n0"bir"1  n« trices,  and  E2,  F2,  F^,  G^,  G^, 

and  H,  are  n  -by-n  matrices.  Since  only  c  and  d  are  needed  to 
x  oo 

calculate  o^g,  we  substitute  (40)  and  (4l)  into  (38)  and  (39)  to 
eliminate  a  and  b.  lhe  two  resulting  relations  can  then  be  written 
in  the  form 


1  -  E2  °1 


c  + 


[-  E2  °J  3  -  [E1  *  E2  <°2  * 


(bk) 


['  F2  0j  3  4  [X  -  F2  G3  -  F3  HlJ  3  =  [Fl  *  F2  <g2  4  g4>  4  f3  H2], 

(^5) 


where  I  is  the  identity  matrix.  If  we  define 


c 

x  =  , 

d 

*  • 

then  (44)  and  (45)  can  be  expressed  as 


(46) 


A  x  =  b  ,  (47) 

and  so  the  desired  solution  is 

x  =  A  ^  b  .  (48) 

We  then  have  the  values  of  c^,  c2,  ...  c^  ,  d^,  d2,  . . .  d  ,  enabling 

o  o 

o  to  be  calculated  by  using  (37)* 

Bb 

All  calculations  were  accomplished  by  using  on  IEM-7090  digital 
computer  programmed  in  FORTRAN.  The  quantities  kb  and  nQ  were  input 
parameters  which  could  be  changed  at  will.  The  sets  of  constants  u, 
v,  B  ,  B^,  and  the  values  of  the  associated  Legendre  funrtions  were 
read  into  the  machine  as  input  data,  whereas  the  values  of  the  radial 
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functions  were  calculated  at  the  beginning  of  the  program,  since  the 

latter  are  dependent  upon  kb.  In  the  case  of  the  iterative  method 

of  solution,  the  values  of  a  ,  b  ,  c  ,  and  d  vere  printed  out,  either 

after  every  iteration  or  after  every  fifth  iteration,  depending  upon 

the  speed  of  convergence  to  the  final  values.  A  subroutine  for  o^g 

was  included  at  the  end  of  the  program,  and  th^  value  of  o  was  also 

Bo 

printed.  In  the  case  of  the  matrix  method,  only  the  values  of  c^, 
d^,  and  were  printed,  as  the  a^  and  b^  were  not  computed  in  this 
latter  method.  Also,  in  the  matrix  method  an  additional  input 
parameter  6q  was  used,  s^  being  the  number  of  terms  retained  in  the 
summations  in  the  elements  of  the  matrices  G^,  G^,  and  H,,  (see 
Appendix  C). 


8.  The  Numerical  Results 

2 

Table  1  lists  60  calculated  values  of  o^g/na  for  50  different 

2 

values  of  ka.  For  seme  values  of  ka,  o__/na  was  calculated  for 

Bo 

several  values  of  n  with  ka  fixed,  In  order  to  determine  the  sensi- 

o  ’ 

2 

tivity  of  Ojg/na  to  n^,  nQ  being  the  number  of  expansion  coefficients 

2 

calculated  in  each  set.  Nine  values  of  a„/na  were  calculated  by 

BS 

using  the  iteration  method.  As  mentioned  previously,  it  was  found 
that  the  iterations  wo.^d  not  converge  for  values  of  ka  above  O.5I8. 


The  matrix  method  was  then  used  for  all  calculations  for  ka  >  O.518. 


Two  different  programs  were  written  using  the  matrix  method.  A 
maximum  of  only  30  coefficients  (nQ  =  30)  could  be  calculated  by  using 
the  first  program,  and  the  maximum  value  of  s^  was  also  limited  to  30. 


22 


TABLE  1 


Calculated  Values  of  the  Normalized  Back-Scattering  Radar  Cross-Section, 


V"8 


2 


0 

Iteration  method 

aB s/nfl2  number  of 

iterations 

Matrix  method 

/  2 

V"a 

e<' 

0 

0.10 

0.0259 

5 

3.362  x  10"6 

50 

0.10 

0.0259 

10 

3.956  x  10"6 

50 

0.10 

0.0259 

15 

4.035  x  10"6 

30 

-6 

0.10 

0.0259 

15 

3.931  x  10 

15 

0.20 

0.0518 

20 

6.395  x  10-5 

50 

0.20 

0.0518 

20 

6.411  x  10-5 

20 

0.50 

0.130 

10 

2.398  x  10'3 

50 

0.50 

0.130 

20 

2.439  x  10_3 

50 

0.50 

0.130 

20 

2.447  x  10" 3 

20 

1.00 

0.259 

25 

3.618  x  10'2 

125 

1.0c 

0.259 

25 

3.614  x  10’2 

25 

1.50 

0.389 

25 

1.591  x  10'1 

25 

2.00 

0.518 

25 

0.4132 

50 

2.00 

0.518 

30 

0.4153 

85 

b  is  the  slant  height  of  the  cone 

(2) 

v  a  is  the  base  radius  of  the  cone 


(3) 

(M 


n  is  the  number  of  expansion  coefficients  calculated  in  each  of 

the  sets  a  ,  b  ,  c  ,  and  d 
v’  u  n  n 


s 

o 


is 


the  number  of  terms  retained  in  the  summations  in  the 
elements  of  the  matrices  Gj^,  and 


o  -> 
L  * 


„ ™ 

0 

Iteration  method 

2 

aBS^nQ  number  of 

iterations 

Matrix  method 

/  2  (4) 

a  rta  s 

BS  0 

2.00 

0.518 

30 

0.4130 

30 

2.25 

0.583 

30 

0.5814 

30 

2.50 

0.648 

30 

0.7716 

30 

2.75 

0.713 

30 

0.9723 

30 

3.00 

0.777 

30 

1 .  lt>+ 

30 

3.30 

0.855 

30 

1-363 

30 

3.60 

0.933 

30 

I.530 

30 

3. 60 

0.933 

33 

1-537 

36 

4.00 

1.04 

30 

1.734 

30 

4.50 

1.17 

30 

2.050 

30 

5.00 

1.30 

30 

2.371 

30 

5.40 

1.40 

30 

2.476 

30 

6.00 

1.55 

30 

2.405 

30 

6.40 

1.66 

30 

2.297 

30 

6.70 

1.74 

30 

2.152 

30 

7.17 

1.86 

45 

1.841 

50 

7.50 

1.94 

45 

1-577 

50 

7.75 

2.01 

45 

1.460 

SO 

8.00 

2.07 

45 

1.412 

50 

8.50 

2.20 

45 

1.462 

50 

9.00 

2.33 

45 

1.495 

50 

9.50 

2.46 

45 

1.367 

so 

10.00 

2.59 

45 

1.077 

50 

10.50 

2.72 

45 

0.8452 

50 

11.00 

2.85 

45 

0.8243 

so 

11.50 

2.98 

45 

0.7746 

50 

12.00 

3-11 

45 

0.5292 

50 

12.25 

3.17 

45 

0.4043 

50 

12.50 

3.24 

45 

0.3273 

so 

13.00 

3.36 

45 

0. 3636 

50 

13.25 

3.43 

45 

0.417b 

50 

13.50 

3.50 

45 

0.4491 

50 

13.75 

3.56 

45 

0.4505 

50 

14.00 

3.62 

45 

0.4490 

50 

14.25 

3-69 

45 

0.5056 

so 

14.50 

3.76 

45 

O.6738 

so 

15.00 

3.88 

45 

1.178 

50 

15.30 

3.96 

45 

1.304 

50 

15.50 

4.01 

45 

1.181 

50 

16.00 

4.14 

45 

0. 70-39 

50 

16. 50 

4.27 

45 

0.8850 

50 

17.00 

4.40 

45 

1.273 

50 

17.  so 

4-53 

45 

0.8129 

so 

18.00 

4.66 

45 

0.7369 

so 

19.00 

4.92 

45 

1.495 

50 

20.00 

5.18 

45 

1.664 

so 
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The  second  program  offered  greater  flexibility,  the  maximum  values  of 

nQ  and  b^lng  45  and  50,  respectively.  These  latter  maximum  values 

of  n^  and  sq  were  limited  by  the  memory  capacity  of  the  IBM-709O  digital 

2 

computer.  Both  of  the  matrix  programs  were  used  to  calculate  o  'na 

Bo 

for  ka  =  0.933,  and  the  values  obtained  were  the  sane,  1.530,  for  the 
same  number  of  terns,  n  .  This  was  reassuring,  considering  the 
completely  dissimiliar  nature  of  the  sequence  of  calculations  in  the 
two  matrix  programs. 

2 

Fig.  3  shows  a  graph  of  a  „/na  vs.  ka,  where  a  is  the  radius  of 

bo 

2 

the  base  of  the  cone  (a  =  b  sin  15°)  •  For  those  values  of  a  /na 

Bo 

calculated  for  several  different  values  of  n  with  ka  fixed,  the  value 
2 

of  Cgg/”0  corresponding  to  the  largest  n^  is  used  in  Fig.  3» 

_  2 

The  values  of  o^/na  in  the  Rayleigh  regior  (ka  <  0.4)  are  not 

2 

shown  in  Fig.  3,  tut  it  is  important  to  note  that  a^/na  obeys 

-4 

extremely  well  the  X  law  predicted  for  this  region.  Furthermore,  by 
using  an  approximate  method  of  Siegel  (Siegel,  1959),  applicable  in  the 
Rayleigh  region,  it  can  be  shown  that  at  ka  =  0.0259  the  normalized 
back-scattering  radar  cross-section  is  approximately  3*56  x  10  The 
value  for  ka  =  0.0259  from  Tfeble  1  is  3  93  x  lO”^*,  vhich  agrees  with 
the  approximate  value  of  Siegel  rather  well. 

For  higher  values  of  ka  the  graph  shows  unexpectedly  rapid 


fluctuations.  It  is  believed  that  these  are  caused  by  convergence 


difficulties,  especially  since  a  curve  obtained  by  using  nQ  =  30  instead 

of  n  =45  showed  even  wilder  fluctuations  for  ka  >  3*2. 
o 


Fig.  3.  Values  of  the  Normalized  Back-Scattering 
Radar  Cross  Section  (c  /no2)  for  n  30-Degree 
Perfectly  Conducting  Cone  of  Slant  Height  b. 
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Also  shown  In  Fig.  3  are  values  of  o_c/n a*"  calculated  by  using 

Bo 

Keller's  modified  geometrical  optics  theory  (Keller,  19oD).  Double 
diffraction  effects  are  incj.uded. 

Mr.  John  E.  Keys  of  the  Defence  Research  Telecommunications 
Establishment,  Ottawa,  Ontario,  very  kindly  supplied  the  present  authors 
with  detailed  data  from  measurements  similar  to  those  on  which  he  and 
R.  I.  Primich  reported  in  the  Canadian  Journal  of  Physics  (  1959). 

Mr.  Keys  has  given  his  permission  for  the  Inclusion  of  these  measure¬ 
ments  in  the  present  report,  and  they  are  plotted  in  Fig.  3-  These 
measurements  were  made  on  flat-based  cones  but  Mr.  Keys  has  informed 
the  present  authors  that  he  has  made  measurements  on  spherically-capped 
cones,  of  the  type  analyzed  in  the  present  work,  and  these  measurements 
are  indistinguishable  from  those  made  on  flat-based  cones. 

Likewise  included  in  Fig.  3  is  the  straight-line  graph  representing 
the  results  obtained  by  using  Siegel's  modified  Rayleigh  theory  (Siegel, 

1*59). 

The  conclusions  to  be  drawn  from  Fig.  3  are  rather  obvious  and  will 
not  be  discussed.  It  is  of  interest,  however,  to  point  out  that  the 
irregularities  appearing  at  ka  values  of  about  1.0,  2.2,  and  2.8  in  the 
curve  illustrating  the  present  work,  occur  in  a  region  where  the  cal¬ 
culated  results  are  believed  to  be  accurate,  so  it  is  considered  that 
these  are  bona  fide  irregularities.  It  was  thought  that  they  might  be 
caused  by  diffraction  from  the  tip  of  the  cone,  but,  when  this  effect 
was  included  i:i  the  Keller- theory  calculations,  the  changes  in  the 
o-values  were  too  minute  to  be  noticeable. 
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In  order  to  look  into  possible  resonance  effects  as  the  cause  of 
these  irregularities,  the  following  table  was  made  up. 


ka 

1.0 

2.2 

2.8 

kb 

3-86 

8.49 

10.82 

2a 

T 

0.318 

0.700 

0 . 892 

b 

\ 

0.6l4 

1.352 

1.723 

na 

T 

0.500 

1.10 

1.40 

b+a 

“T 

0.773 

1.702 

2.17 

The  values  of 

na/\  make  it  appear 

that  these 

irregularities  may  be 

resonances  in 

response  caused  by  current  paths  1‘iom  top  to  bottom  of 

the  cone  along  the  edge  of  the  base.  In  viev  of  the  fact  that  the 
base  edge  of  the  cone  is  very  important  in  determining  the  scattering 
characteristics,  it  is  not  surprising  that  these  resonance  effects 
occur. 
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APPENDIX  A 

ANALYTIC  SOLUTION  FOR  EXPANSION  COEFFICIENTS 


The  solution  for  the  expansion  coefficients  will  now  be  presented 
in  its  analytic  detail. 

1  dpl 

First  we  multiply  (27)  by  P^(cos  0) ,  multiply  (28)  by  sin  9  , 

m  do 

and  subtract  the  two  results,  obtaining, 


l  Vn<kb) 

n*l 


r1  p1 


m  n  +  sin  9 


dp1  dP1 


sin  9 


+  dP*  I 
m  n 

- 


d9  d0 


"  r  dP1  dr1 

♦  )  r  J '(kb)  I  P1  ^  +  P1  1 

La  n°n'  _  m  d0  d0  n  J 

n=l 

"  r  ^  ^  dP1  dP1  , 

+  )  c  h  (kb)  ‘  ~ |  +  sin  9  -r~ 

L  n  n'  L  sin  0  d9  d0  J 

n=l 

®  dP1  dP1 

♦  )  d  h'(kb)  P1  TT-  +  T32  P1 

La  n  n  L  m  d9  d0  n  J 


n=l 


r 


)  °v>v<kb> 

"  .J  ^  A  J  A  JA  “J 


sin  0 


dd  d0 


< 


l  bu-c<kb)  [  *i  Si  *  fa  ].  °  < e  <  ® 


de  de 


0,  9  <  9  <  it  . 

’  o  - 


(1-1) 


By  integrating  both  sides  of  (l-l)  with  respect  to  9  over  the  interval 
0  to  n  and  combining  terms,  we  obtain, 


V  P  1  "  r  P1  l'1 

)  y  J  (kb)  >  c  h  (kb)  I  m  n 


n  n  n  n  J 


n=l 


0 


sin  0 


.  n  dP1  IF1 
sin  0  m  n 

dO~  do“  J 


dO 


)  r  r  j'(kt.)  ,  d  h'(kb)  1  j  ”r  p1  <  ,  K  p1 

L.  ^  n  n  nr.  j  J  m  t-  tt  n 


dO 


n=l 


dO  dO 


'  l°vJv(kb)  J 


0  . 
r  °r  p1  p 


m  v  +  sin  0 


dF* 


^  dF7  1 

n  v 


dO 


sin  0 


dO  dO 


♦  )  bj'(kb)j  T 

L  u°u  J  L  m  tt  t—  u  J 


dO  . 


(1-2) 


0 


dO  dO 


The  first  integral  in  (1-2)  is  n  familiar  integral  of  boundary  value 
problems . 


P1  P1  .  J  „  dP1  dP1  1  _  2  m(m*l) 


0 


m  n  sin  0  _ ra  _ n  dO  = 

sin  0  dO  dO 


T 

J  6 


2nr+l 


m  n 


(1-3) 


where  5  is  the  Kronecker  delta.  The  second  and  fourth  integrals  ore 

m  n  B 


easily  evaluated  by  using 


P*  I  -  0  ,  P1  | 


n  .  -  0  ,  p;  !  =  o  ,  and  y1  I 

n  0-0  0-n  U  0=0 


r  r  P1  .  <  P1  1  dO  .  I'  d  V  P1T  =  Tpl  pll  |”  .  o 

J  L  m  37—  377-  n  J  J  ..  m  n^  ^  m  nJ' 

0  dO  dO  Q  0 


(1-M 


r  °r  pi  dp^  +  dp^  pi  ■  d0  =  drpi  pi]  =  rpi  pi-!  j  c 

J  .  m  375  375—  ^  j  J  -  m  uJ  .  m  uJ 

0  ie  do  0 


=  0 


dO  dO 

(J  V. 

(1-5) 

The  third  integral  is  eva]  uoted  by  utilizing  the  associated  Legendre 


equation, 


(1-6) 


d  /  \  11 

—  sin  0  n  t  +  n(rH-l)  sin  8  F  -  P 
dO  V  -r—  n  n 

d0  - —  ^ 

sin  0 


Multiplying  by  P1  nnd  using  the  product  differentiation  rule,  we  obtain 

m 


P1  F1 
8 


dp1  dP1 


5 — 8  *  sia  0  _ E  _ —  -  5pn  q  ^  ..  n(n+l)sin  fl 

in  9  d 0  dO  dO  V.  m  -r^  J  n^l^sln  9  pm  \  • 

dO 

(1-7) 


Integrating  free  0  to  9  ,  there  results 


9 


I  ♦  Sir.  8  K  51d0  .  Fein  e  P1  51  r-  r  8  P1  P1 

0  "sin  6  ,n  J  L  n  J  ~  m  n 


dP1  dp1' 
_ m  _ r 

dO  dO 


,  e  e 

dpl->  o  o  o 


dO 


0 


0 


(1-8) 


The  integral  in  (1-8)  con  easily  be  evaluated  by  using  (l-7).  If  m  and 
n  are  interchanged  in  ( 1- 7 )  and  the  result  subtracted  from  (1-7),  the 
ensuing  equation  is 

rm(nr+l)-n(n+l)  sin  0  P1  P^"  =  ~  ^  sin  0  P^  ^n  -  sin  0  P1  ^nl  . 

m  n  dO  _  m  -  "  - 1 


The  integral  in  (1-8)  may  now  be  evaluated. 
0 


m 


dO 


n  de  J 


(1-9) 


J  n(  n+l)  sin  0  P1  P1  dO  =  t  rsin  9  (p1  ^  -  F1  dP^)"l| 

JQ  =  r‘  m(m*l)-n(n*l)  -  o  V  m  n  Q 

(1-10) 


Substituting  (1-10)  in  (1-8)  and  combining  terms,  there  results 


The  third  integral  in  (l-2)  is  then  obtained  by  letting  n  =  v  in  (l-ll) 
and  making  use  of  the  condition  (26). 


de. 


«,  < 


'  '  r  P1  ^  A  .  .  dP1  dp1  .  -  v(  v+ 1 )  sin  0  "1  dr1!, 

m  v  +  sin  0  _ m  _ v  idO  =  _ o  P  m  ( ! 

* n  sin  0  dO  dO  J  m(  m  1 )  -  v(  v+l)  v  dO  “ 


(1-12) 


Upon  substituting  the  integral  resalts,  (1-3),  (l-u),  (1-5),  and 
(1-12),  into  (l-2),  we  obtain 


r  "  *- 

)  fy  J  (kb)  4  c  h  (kb)"l  2  m(m+1)j  5 

L.  l  n  n  n  n  J  - ; - - - 


2m  1 


v(vtl)  sin  9  r  ,  1-. 

*  -1  Jji  o  > 

dO  o 


(1-13) 


which  can  be  solved  for  c  : 

m 


c  =  -  Y  J  (kb)  *  (2ml)  sin  0 

m  mm  o 


dP1 

— El 


Vkb)  2  m(ml)  2  h  (kb)  d°  ^°c 

j  m 


%  v('"l>Jv(kb)I4(  cos  0^) 
v(  nh  l ) -m(  m  1 ) 

(l-l1*) 


The  separation  of  d^  is  accomplished  in  a  similar  manner  by 

dpd  1 

multiplying  (27)  by  sin  0  n  ,  multiplying  (28)  by  p  ,  and  subtracting 

dO"  '  m 

the  two  results,  obtaining, 


w 

)  V  j  (kt)  r  r  ^  <  P1  1 
4-  n  nx  _  ra  -r—  n  J 


dO  dO 


CD 

\  r  4  0  dP1  d^1  1 

+  )  1  J  (kb)  m  n  +  sin  0  m  n  1 

,  n  n  L  sin  0  dO  dO  J 

n=l 

CD 

)  c  h  (kb)  r  K  F1  .  r  iFn  1 

L.  n  n  _  -rr—  n  m  -r—  j 


\  ,,  l'/., \  r  p1  p1  .  „  dp"  ip^  1 

)  d  h  (kb)  ra  n  ■*  sin  y  m  n 
l-i  n  n  -  ■  . _  r  •>  ^  —  J 


dP1  dP1 


sin  0 


dO  dO 


3^ 


r 


=  j 


i  .,„<»>  t  <  <•  <><  ] 

V 

4 1  Vi(kb)  r  i  *  Sin  S  ^  ]  -  0  ;  8  <  «0 

U  U  J.  L  nir,  a  1,1,1  do  d© 


0  ,  0  <  ©  <  n 

o  —  — 


(1-15) 


Integrating  both  sides  of  (1-15)  with  respect  to  0  over  the  interval 
0  to  n  and  combining  terns,  we  obtain, 

)  y  J  (kb)  ♦  c  h  (kb)]  I  r  P1  dFt  +  ^  P1  ]d© 

_  n  n  n  n  J  J  m  -rr-  -rr—  n  J 


n=l 


d©  d© 


“  n 

♦  )  r  T  J^kb)  +  d  h^kb)^  ^  F  ^n  +  sin  0  ^n  d© 
t->  {.  n  nx  n  n  'J  JA  L  -nr-  -rr-  J 


n=l 


0  sin  0 


d©  d© 


*1  Vv<kb>  J + 


+  ^  b^jJ(kb)  j  I"  ' m  ^  ■*■  sin  0  ]d©  .  (l-l6) 


0  sin  6 


d  ©  d© 


The  values  of  the  first  and  second  integrals  in  (l-l6)  are  given  by 
(l-U)  and  ( 1- 3) >  respt  ctively.  The  third  integral  is  easily  evaluated. 

©  ,  e 


J  r  dPm  pj;  ♦  p1  dPv  ]d©  =  f  0  d,  P1  p1]  =  p1  (cos  ©  )  pj;  (cos  0  ) 
J  L  -T7T-  v  m  tt  j  L  m  vj  m  '  o'  v  '  o' 

(1-17) 


d© 


de 


'0 


The  fourth  integral  is  evaluated  by  using  (l-ll)  with  n  replaced  by  u. 

Q 

r  0  r  P1  P1  .  ,  „  dP1  dP1  m(m+l)  ,  „  frX  dP1!  1 

m  u  ♦  sin  0  m  u  d©  =  -j — =4 — T  -*  sin  0  r  u 

Jo  l7uTe  !TW-j  ouo^JI^ 

(1-18) 
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By  substituting  these  results  into  (l-l6),  we  obtain 
00 

J  .  1? 

I  T  j  ( kb )  +  d  h~(  kb)  I  ?  ' 
l  n  n  n  n  J  - 

n=l 


)  rrj:(kb)  ♦  d_h^kb)l  2[m(m4l)]2  5 

L,  l  n  n  n  n  J  — : -  mn 


2m  1 


*  1  %Jv(kb)  pm(c°3  “o'  Pi(cos  “o' 


m(ml)  sin  Q  r  1  n 

)  b  j'Ub)  - yr - j - -y  rp1  dPU  ll 

L  u°u  m(mlj-u(u+l)  L  m  J  e=e  ' 


(1-19) 


which  can  be  solved  for  d  : 

m 


d  = 

m 


r  j'(kb)  (2ml)  (cos  9  ) 
■  I_  m  +  _ m _ o 

2|m(mT)ph^(kb) 

j  m 


m  r~f 


hTkb T 

IQ 


)  av^v(kb^  ^  (coS  6 J 


(2ml)  sin  G  PP  (cos  6  )  -  b  i'(kb) 


o  m 


2m( ml )  h'(kb) 


m 


2-  )  - 
L>  m 


dP 


(mi)-u(u+i)  de 


d=Q 


(1-20) 


The  separatlrn  of  is  accomplished  by  multiplying  (29)  by 


dP  i 

sin  9  _ a  >  multiplying  (  30)  by  and  adding  the  results,  obtaining, 

dO 


O' 


>  Ynj'(kb)  [  ^  ♦  sin  0<  K  ] 
L*  nn  b  _ .  _  ~  j-  0 


n=l 


sin  9 


d0  d0 


.)  r  j  (kb)  r  p1  < ♦  «£  p1] 

L  n  n  L  a  -rr-  n  J 


n=l 


de  de 


00 

*  )■  c  h'(kb)  r  ^  *  sin  0  «£  1 

L  n  nv  L  rr— ^  -r sr-  ttt-  J 


n=l 


sin  e 


de  de 


J6 


oo 

-  )  d  h  (kb)  T  p1  ^  ^  p1  1 

L  n  n  L  a  -rr-  rr  n  J 


n=l 


d  0  dd 


)  .!>>  [  £  ^  *  sin  9  <  ] 

V  v  L  _  4  _  r\  jTs\  j  a  -J 


sin  0 


dd  dd 


-)  b  j  (kb)  [  P1  ^  +  ^P1  ]  ,  0< 
L  uV  L  or  ITT*  TS“  u  J 


dd  dd 


d  <  d  . 
o 


(1-21) 


where  the  subscript  or  denotes  a  particular  value  of  the  Infinite  set 
of  v.  Integrating  both  sides  of  (l-2l)  with  respect  to  d  over  the 
interval  0  to  dQ  and  combining  terms,  we  obtain, 


1  [  Vn(kb)  +  cnhn(kb)]  I  °  [  -2-Ia  +  flln  »  ^  ^  ]  d» 


n=l 


sin  d 


dd  dd 


-  I  [r„J„<kb>  *  dnhn<kb>]  I  °  K  ^  ^  ]  d8 


n=l 


dd  dd 


- 1  %jv(kb>  j  °  [  + sin  a  ^  ^  ] 


dd 


sin  d 


dd  dd 


iu'“>  ■ 


(1-22) 


dd  dd 


The  first  integral  in  (1-22)  is  evaluated  by  using  (l-ll),  replacing 
m  by  cr  and  noting  that  or  is  a  particular  value  of  the  set  of  v. 

Q 

r  ,  .  dp1  ar1 1 or(of^l)  sin  d  r^-ii 

or  n  ♦  sin  d  or  n  dd  *  '  o  r  n  | 

J  L  sin  e  ~ST  d8  J  a(a*l)-n(n+l)  L  0  dS  J  0*9° 


(1-23) 
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The  second  and  fourth  integrals  are  easily  evaluated. 


J  Fp1  ^  +  ^  p1  ]  d0  =  |  df P1  P1 1  =  P^(cos  0  )  P^(cos  0  ) 

J  L  or  -rr—  S7 r"  n  J  j  I.  a  nJ  or  o  n  o 


d0  d0 


(1-24) 


6  & 

|  °  Tp1  <  *  p1  1  ae .  [  °  ^p1  p1]  .  iV  r1]! 

J  L  or  -r?r-  -T77-  3  l.  or  uJ  La  uJ1 


d0  d0 


(1-25) 


The  third  integral  is  evaluated  by  using  (l-ll), 


0  i  i 

r  °r  r  r 

j  ^  at  v  +  sin  0 


sin  0 


dP1  dP1  I 
O'  _ v 

d0  d0  J 


d0  e  5 


F  +  ain  0  (  dJoT)  1  d0  , 


6in  0 


(1-26) 


where  6^v  is  the  Kronecker  delta.  The  integral  in  (1-26)  can  be 
evaluated  by  the  use  of  (1-7).  IT  (1-7)  Is  re-written  with  m  and  n 
replaced  by  a,  there  results 


+  sin  0  ^  ~  (^sin  0  P^  ^  +  o(o^l)  sin  0  (P^) 


sin  0 


(1-27) 


ht  integral  in  (1-26)  is  then 


J  °  F  +  sin  0  ^  ^er)  ]  d0  =  (sin  0 


sin  0 


d0  o 


g(on-l)  I  sin  0  (P1)  d0  .  (1-28) 

j  at 


The  first  term  on  the  right  side  of  ( 1-28)  is  seen  to  vanish,  and  then, 
if  we  define 


B  =  [  sin  0  (P1)  d0  , 
a  J  at 


(1-29) 


ifi 


for  (1-26)  there  results 
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r  o  r  y~^^~  j ^.1  ..A  -* 

I  or  v  +  sin  0  cr  v*  d  9  =  or(o^l)  B  6  .  (1-30) 

-  u  _  j  _  n  jrt"  jrt  -J  or  O' v 


sin  d 


dd  dd 


Upon  substituting  these  results  in  (1-22),  we  obtain 
00 

5  r  v  rub)  *  =  h'(kb)i <,(<r,i)  sin  ®o  r  p1  dpn  1 1 

L,  1  "  ”  "  n  J  or(a*l)-n(n-»l)  L  “  W  J  B*B 

n*l  c 


Hr  J  (kb)  >  d  h  (kb)l  P1  (cos  9  )  P1  (cos  0  ) 
n  n  n  n  J  or  '  o  n  o 


n=l 


1 


a  Jy(kb)  Q'(o^-l)  B  6  , 

v°v'  '  '  7  o  crv  ’ 


(1-31) 


and  solving  for  a^  , 


w 

oln  9„  pi  (cob  9^  <  Ynj'(kb)  -  c.hl(kb)  d/ 


o  cr 


n  n 


Bcr  n=i  -  n(n+l)  dd 


d=d 


P1  ( C08  0  ) 


cr(or^l)  B  y 


Q  ) 

rjW  I  [rnJn<kb)  *  dnhn(kb)]  *£  <cos  8o>  • 
O'  cr  , 

n=l 


(1-32) 


The  last  coefficient  to  be  separated  is  b^.  This  is  accanplished 


dd 


.  dP1 

by  multiplying  (29)  by  P^  ,  multiplying  (30)  by  sin  d  ft  ,  and  adding 

the  results,  obtaining, 

00 

dr1  .  dPi 


Lv>) 

-[  rnJn(kb)  [  *  Bln  9  5  5  ] 


n=l 


sin  d 


dd  dd 


X  cnhn(kb)  r  ^  ^  +  ^  P1  ^ 

u  n  n  -3  T7T  ITT  n  J 


-  )  d  h  (kb) 
L  n  n 


H  dO  dO 

P  P*  P1  .  .  .  dph  dP1 
13  n  +  sin  0  ft  n 

sin  0  d0  d0 


■  i  %x<»>  [ -i  £  •  3 ’i  ] 


d0  dO 


-  )  b  J  (kb)  T  Pft  -►  sin  0  ^  ^  ^  ,  0  <  0  <  0  ,  (l 

L*  u  u  L  ■  11  t  ~rr~  tt  J  -  o 


■  in  0 


dO  dO 


where  the  subscript  ft  denotes  a  particular  value  of  the  infinite  set 
u.  Upon  integrating  both  sides  of  (1-33)  vith  respect  to  0  over  the 
interval  0  to  0^  and  combining  terms,  we  obtain, 

)  Y  ♦  c  h'(kb)l  f  r  pj;  ^  +  dFp  P1  1  d9 
L.  L  "  '  n  n  -  J„  -  e  W  iO  "  J 


1  [  fnJn(kb)  *  dnhn(kb)  ]  J  [  ♦  sin  8  "g  “n  ]  de 


sin  0 


dO  d0 


■l  v;<«»>J  K^-3k]» 


d0  d0 


•I  bA(kb)je0[i^.si„9^<]de.  (1.3b) 


sin  0 


dO  d0 


The  first  integral  in  (1-34)  is  easily  evaluated,  noting  that  3  is  a 
particular  value  of  the  set  u. 


(1-35) 


The  second  integral  is  evaluated  by  using  (l-ll)  and  replacing  m  by 
3,  obtaining, 


(1-36) 


The  third  integral  in  (1-3*0  also  is  elementary. 


(1-37) 

The  fourth  integral  is  evaluated  by  using  ( I-  .  1 ) , 


(1-38) 


where  5^  is  the  Kronecker  delta.  Phe  integral  in  (1-38)  is  evaluated 
by  using  (1-28)  with  a  replaced  by  (3. 


Th«  first  term  on  the  right  side  of  ( 1- 39 )  is  seen  to  vanish,  and  then, 
if  we  define 

*o 

Bp  =|  sin  0  (f£)  d8  , 

0 


(1-40) 


for  (1-38)  there  results 


r  p1  p1 

L  -§-±L  + 

sin  0 


sin  0  "Iu0  =  0(^+1) 

r\  j  r\ 


dO  d0 


B0  &0^ 


(1--1) 


Upon  substituting  the  results  of  these  integrations  into  (1-34).  we 
obtain 


■l  [riA(lrt> 


n=l 


+  d  h_(kb)l  n(nU)  sln  So 

“  n  J  n(n+l)-0(0r+l) 


•■l  H(kb>  Be V 


(1-42) 


and  solving  for  bQ  , 

0 


sin  0 

o 

’p  '  e(p*D  BpjBUb) 


dp1 

Zfi  I  ) 

de  'M  L, 

o  n=l 


[rnJnUb)  ♦  dnhn(kb)] 


n(  n+l) -0(0+1) 


n(  n+l)  P^( 
(1-^3) 


Equations  (l-l4),  (l-20),  (I-32),  and  (1-43)  represent  the  formal 
solution  for  the  expansion  co-efficients,  and,  when  the  values  of 
and  r*n  are  substituted  by  using  (8),  are  equivalent  to  (31)  through 
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APPENDIX  B 

LEGENDRE  FUNCTION  CONSTANTS 


u(l)(2)  v(l)(3)  b(M  b  (M 

u  v 


1 

1.031631 

0.967140 

1.310 

1.35806 

2 

2.084434 

1.918899 

2.346 

2.42491 

3 

3.1^9929 

2.887078 

3.3^7 

3*379^5 

4 

4.223096 

3.887853 

4.341 

4.28564 

5 

5 . 301087 

4.917100 

5.3^0 

5.18033 

6 

6.382249 

5.965629 

6.323 

6.09038 

7 

7.465580 

7.026428 

7-168 

7.03236 

8 

8.550454 

8.095125 

8.206 

3.058 

9 

9.63645 

9.169073 

9.227 

9.009 

10 

IO.72329 

10.24665 

10.24 

9.973 

11 

11.81078 

11.32681 

II.25 

10.95 

12 

12 . 89879 

12.40890 

12.25 

11-93 

13 

13.98718 

13.49242 

13.26 

12.91 

14 

15.07592 

14.57706 

14 . 26 

13.89 

15 

16.16491 

15.66258 

15.26 

14.88 

16 

17-25414 

16.74882 

16.26 

15.87 

17 

18. 3^35^ 

17.83562 

17.26 

16.86 

18 

19- ^3311 

18.92291 

18.26 

17.86 

19 

20.52280 

20.01059 

19.26 

18.85 

20 

21.61262 

21.09860 

20.26 

19-84 

21 

22.70252 

22.18690 

21.26 

20.84 

22 

23.79253 

23.275^5 

22.26 

21.83 

23 

24  88260 

24.36421 

23.26 

22.83 

24 

25.97275 

25.^5315 

24.26 

23.83 

25 

27.06294 

26.54226 

25.26 

24.82 

26 

28.15320 

27.63151 

26.26 

25.82 

27 

29.24349 

28.72088 

27.26 

26.82 

28 

30.33385 

29.81037 

28.26 

27.82 

29 

31.42421 

30.89995 

29.26 

28.81 

30 

32.51465 

31.98963 

30.26 

29.81 

(1) 

Donated  by  Dr.  P.  C. 

Waterman 

of  AVCO. 

(2) 

Determined  from  P^( 

cos  165°) 

=  0  . 

^  ^  Determined  from  dP^  =  0  . 

dT^e=i65° 

(h)  r  1650  1  2 

'  '  Defined  by  =  J  sin  0  (P^)  d0  . 


0 
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31 

33.61 

33.08 

31.26 

30.81 

32 

34.70 

34.17 

32.26 

31.81 

33 

35.79 

35.26 

33-26 

32.804 

34 

36.88 

36.35 

34.26 

33.803 

35 

37-97 

37.44 

35-26 

34.801 

36 

39.00 

38.53 

36.26 

35.800 

37 

40.15 

39-62 

37.26 

36.798 

30 

41.24 

40.71 

38.26 

37.797 

39 

42.33 

41.80 

39.26 

38.796 

40 

43.42 

42.89 

40.26 

39.794 

41 

44.51 

43.98 

41.26 

40.793 

42 

45.60 

45.07 

42.26 

41.792 

43 

46.69 

46.16 

43.26 

42.791 

44 

47.78 

47.25 

44.26 

43.790 

45 

48.87 

48.34 

45.26 

44.789 

46 

49.96 

49.43 

46.26 

45.788 

47 

51.05 

50.52 

47.26 

46.787 

48 

52.14 

51.61 

48.26 

47.787 

49 

53-24 

52.70 

49.26 

48.786 

50 

54.33 

53-79 

50.26 

49.785 

51 

55-42 

54.88 

51.26 

50.784 

44 


n 

P^(cOS  165°) 

F^(cos  165°) 

dp1 

ni 

dO  1 0*165° 

dP1 

4  | 

dO  1  0=165° 

1 

-0.25881924 

-0.52346792 

0.96592579 

1.88b 

2 

0.750 

1.4156256 

-2.5980742 

-5.177 

3 

-1.4228831 

-2.3050079 

4.3396881 

9.503 

4 

2.2069309 

2.98087 

-5.4575812 

-14.65 

5 

-3.0177961 

-3.50544 

5-1518426 

20.48 

6 

3-7646396 

3.92272 

-2.6831313 

-26.91 

7 

-4.3581639 

-4.27745 

-2.556614 

32.93 

8 

4.748 

4.603 

11.19 

-40.63 

9 

-4.797 

-4.891 

-22. 08 

48.73 

10 

4.514 

5.159 

35-26 

-57.23 

11 

-3-884 

-5-413 

-49.72 

66.12 

12 

2.922 

5.655 

64.10 

-75*39 

13 

-1.675 

-5.880 

-76.74 

85.03 

14 

0.2149 

6.109 

85.84 

-95.03 

15 

1.364 

-6.323 

-89.66 

105.4 

16 

-2.953 

6.530 

86.71 

-116.1 

17 

4.433 

-6.732 

-75.92 

127.I 

18 

-5.690 

6.927 

56.81 

-138.4 

19 

6.618 

-7.117 

-29.64 

158.0 

20 

-7.133 

7.302 

-4.551 

-162.0 

21 

7.176 

-7.483 

43.92 

174.2 

22 

-6.720 

7.659 

-85.90 

-186.7 

23 

5.775 

-7.631 

127.3 

199.5 

24 

-4.388 

8.000 

-164.8 

-212.6 

25 

2.637 

-8.166 

194.7 

226.0 

26 

-0.6337 

8.328 

-213.6 

-239.6 

27 

-1.491 

-8.487 

218.7 

253.5 

28 

3-591 

8.643 

-208.2 

-267.6 

29 

-5.516 

-8.796 

180.8 

282.0 

30 

7.126 

8.947 

-137-1 

-296.6 

31 

-8.296 

-9.096 

78.73 

311.5 

32 

8.929 

9.242 

-8.613 

-326.6 

33 

-8.964 

-9.385 

-68.98 

342.0 

34 

8.380 

9.527 

148.8 

-357*5 

35 

-7-199 

-9-666 

-225.2 

373.4 

36 

5-487 

9.804 

291.9 

-389.4 

37 

-3.349 

-9-940 

-343.2 

405.7 

38 

0.9210 

10.07 

374.0 

-422.1 

39 

1.634 

-10.21 

-380.2 

438.8 

4o 

-4.142 

10.34 

359.5 

-455.7 

4l 

6.427 

-10.46 

-311.2 

472.9 

42 

-8.324 

10.59 

237.0 

-490.2 

43 

9.693 

-10.72 

-140.3 

507.7 

44 

-10.42 

10.84 

26.61 

-525.5 

45 

10.45 

-IO.96 

96.96 

543.4 

4o 

-9.765 

11.09 

-222.1 

-561.5 

47 

8.388 

-11.21 

339.7 

579.9 

48 

-6.404 

11-33 

-440.9 

-598.4 

49 

3-937 

-11.44 

517.1 

617.1 

50 

-1.149 

11.56 

-561.4 

-630.O 

51 

-1.773 

-11.68 

568.5 

655.1 

APPENDIX  C 


DEFINITIONS  OF  MATF.IX  ELEMENTS 


The  elements  cf  t:.e  me  trices  ir  (36)  through  (4l)  ere  given  by 
(31)  through  (  3M  •  If  we  let  [A]  deriOte  the  J t'f ‘  eiement  of  the  1^ 
row  of  a  matrix  A,  then  the  matrix  elements  ere  defined  by 


M 


ir‘  (2m*  1)  Jr(kb) 
m(m*l)  h  (kb) 


(3-’) 


m 

L  2. 


(2ml)  sin  eoJ^\0__Q  v  ;VD  *v  (kb)  Pv  (cos  0  ) 

_ c _ a _ r _ 

2^"rr.(  rr*  1  ) ^  h  ( kb  )  v  (  v  1 1 )  -  m(  nr*  1 )  ”1 
L  J  m  ^  n  n  J 


(  3-D) 


M 


i^1  (2m*l)  J^(kt  ) 

m(m+l)  h'(kb) 
m 


(3-3) 


(2m*l)  (cos  o  )  j  (kl)  P*  (cos  0  ) 


2  n.(  rr*  1 )  |  n  '  (  kb ) 
u  J  rr. 


(3-M 


(2m*l)  sin  6o  P^  (  ccs  0^)  T  (kb) 


o'  “u 


do  'o-e 


2r,(  m*  1 )  r. '(kb)  ^ rr*  1 )  -  u  (u  +i)l 
rr.  w  r.  n  j 


(3-5) 


rG  1 

tji. 


rG  1 

w  d  _ 


sir.  »c  Pv  (cos  6o)  h^Kb)  35-^ 


JC  (’fth)  *v  V1)  -  n(ml) 

rr,  r. 

S  dp1 

sir.  0o  r‘  (cos  Oo)  r  ln  (2„U)  J^(Kb) 


'  '  (kb  T) 

*•  v  KC  -  v 


,  n(  n*  l)  v  (  v  +  1 )  -  n(  n+  I ) 

n-  j.  u  rr.  m 


(3-6) 


(3-7) 


(3-8) 


-  pj;  (cos  6  )  h  (kb)  P*  (coe  6  ) 
v  on  n  o 

m 

vn<v1)Bv 

m  m 


p1,  (cos  e) 


m 


l 
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n  n 
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m  m 


dPJ 


n(  n+l) 


(3-9) 


oin  eo  d 0~'e»e  hn(kb)  n(n4l)  Pr.  (c°8  0O} 

_ o 
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n  m 
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U  and  v  represent  the  value  of  the  sets  u  and  v,  respectively, 
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